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Objectives  of  Research  Proornm 

"“^The  research  aimed  at  fir.dinp  and  developing  methods  for  estimating  the 
main  features  of  response  of  enpineerinp  structures  subjected  to  severe  dynamic 
loading  of  pulse  type,  with  emohasis  put  on  methods  valid  both  for  lame  deflec- 
tions and  for  structures  of  materials  exhibitinp  stronp  strain  rate  sensitivity 
in  the  plastic  ranpe.  Problem  types  of  practical  importance  include  explosive 
loading,  either  external  due  to  military  attack  or  internal  for  example  due  to 
disruptive  accident  in  a pressure  vessel  or  containment  structure;  various  tvpes 
of  vehicular  impact;  wave  impact  on  ship  or  offshore  structures;  and  hiph  enerpv 
rate  forminp.  Preliminary  applications  have  been  made  of  the  methods  under  in- 
vestipation  in  the  propran  presently  beinp  reviewed  in  all  but  the  last  of  the 
above  areas. 

The  ne^c^for  reliable  estimation  tochninues  arises  in  part  simnlv  from  the 

\ 

need  to  extend 'basic  knowledpe,  and  in  part  from  the  inadeouaev  of  approaches  by 
wholly  numerical  schemes.  The  latter  have  been  the  object  of  research  and  de- 
velopment on  a large  scale  over  the  past  decade,  particularly  for  methods  of 
finite  element  tvpc.  Desnito  this  effort  examples  are  familiar  to  most  of  us 
of  cases  where  a supposedly  tested  numerical  code,  anpliod  in  slipht.lv  different 
circumstances,  lias  led  to  misleadinp  or  totally  nonsensical  results,  often  after 
substantial  computinp  cost.  It  is  not  surprisinp  that  this  is  the  case.  The 
problems  concerned  are  intrinsically  difficult  when  approached  hv  numerical 
methods.  These  are  necessarily  of  stop-bv-step  nature,  proccodinp  from  pivon 
initial  surface  tractions  and  continuing  with  specified  loadinp  and  constraint 
conditions.  The  initial  response  consists  of  propagation  of  various  types  of 
elastic  wavss.  Traverses  of  these  waves  lead  to  plastic  strain  waves,  which 
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build  up  large  deformations  from  a complex  interplay  of  strain  rate  dependent  and 
irreversible  dissipative  mechanisms.  Problems  of  convergence  and  stability  are 
a fundamental  hazard  in  all  numerical  schemes,  which  can  probablv  never  be 
eliminated.  The  need  for  simple  but  reliable  and  rigorously  based  estimation 
techniques  is  obvious,  both  as  an  adjunct  to  full  numerical  methods,  and  as  a 
means  of  gaining  better  understanding  of  basic  response  phenomena  that  wholly 
numerical  methods  are  incapable  of  providing. 

The  estimation  methods  investigated  under  the  program  here  reviewed  were  of 
two  types:  theorems  giving  bounds  on  deflections  and  response  times,  and  methods 
in  which  "natural''  responses  in  simple  mode  (separated  variable)  form  are  used 
to  approximate  the  actual  dynamic  response  history.  These  are  based  directly  on 
overall  energy  and  momentum  conservation  principles:  the  bound  theorems  on  the 
minimum  potential  energy  theorem,  and  the  mode  approximation  technique  on  a 
"minimum  error"  device  which  can  be  identified  with  a statement  of  momentum  con- 
servation. Such  approaches  can  be  expected  to  guide  and  strengthen  the  analyst's 
intuitive  feel  for  the  essential  features  of  the  response.  While  the  numbers  ob- 
tained will  not  be  exact  or  in  some  cases  even  good  approximations,  thev  will  be 
cheap  and  will  rarely  be  misleading  or  nonsensical,  which  cannot  be  said  for 
wholly  numerical  approaches.  Although  the  problems  dealt  with  so  far  are  proto- 
type problems  of  fairlv  modest  complexity,  results  from  them  generally  indicate 
that  the  objectives  stated  are  realistic,  even  though  thev  have  not  vet  been  fullv 
attained. 
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Work  Done  Including  Reports  and  Publications 

The  results  obtained  are  described  in  nine  reports.  These  are  the  first 
nine  items  in  the  list  of  Reoorts  .and  Publications,  and  referred  to  here  by 
number  in  Part  A of  that  list.  Some  additional  references  cited  are  included 
in  Part  3 of  the  same  list. 

The  two  approaches  investigated  referred  to  as  Deflection  Bounds  and  the 
Mode  Approximation  Technique,  although  conceptually  distinct,  are  in  fact  linked 
to  each  other.  However,  for  convenience  work  is  summmarized  under  the  two 
headings . 

Deflection  Pound 

A new  bound  theorem  was  presented  Cl,  2]  for  a structure  subjected  to  im- 
pulsive loading.  The  theorem  is  applicable  to  structures  whose  material  exhibits 
strongly  rate  sensitive  behavior  in  the  plastic  range,  and  applies  to  finite 
strains  and  finite  deflections.  The  latter  are  of  particular  importance  when 

I 

the  structure  geometry  and  constraints  are  such  that  substantial  qualitative 
changes  occur  in  the  stress  field  during  the  response  and  the  new  theorem  allows 
consideration  of  these  effects  together  with  viscoplastic  behavior. 

The  new  theorem  is  based  on  the  principle  of  minimum  potential  energy.  We 
are  considering  a structure  with  specified  kinematic  constraints  ,\nd  initial 
velocity  field  that  corresponds  to  the  applied  impulsive  pressures.  The  material 
behavior  mav  be  of  anv  type  including  time  dependent  inelastic  as  well  as  elas- 
tic behavior.  Tor  such  a material  (nonholonomic)  the  potential  energy  minimum 
principle  is  not  valid,  in  general.  It  becomes  valid  and  is  available  in 
general  form  for  our  purposes  when  use  is  made  of  the  concepts  of  "extremal  paths" 
in  strain  and  stress  space  which  render  minimum  and  maximum  the  specific  work 
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and  specific  complementary  work,  respectively.  These  concepts  provide  the 
rigorous  basis  of  the  new  theorem.  A static  load  traction  svstem  is  introduced. 
The  bound  on  final  deflection  at  anv  point,  due  to  the  dynamic  loading,  is  ob- 
tained in  terms  of  the  deflections  reached  under  the  defined  static  loads.  The 
total  strain  energy  in  the  structure  due  to  the  static  loading,  computed  ac- 
cording to  extremal  path  material  behavior,  must  be  at  least  as  great  as  the 
initial  kinetic  energy  corresponding  to  the  specified  initial  velocity  field. 

In  applying  the  theorem  to  several  prototype  problems,  as  in  [1,  2,  8],  no 
attempt  was  made  to  use  constitutive  laws  expressing  general  time  dependent  be- 
havior. (These  are  not  well  established,  in  anv  case,  since  the  required  ex- 
perimental results  are  not  available).  The  simplest  forms  appropriate  for  the 
more  important  stronglv  rate  sensitive  metals  such  as  mild  steel,  stainless 
steel,  and  titanium,  are  ones  of  essentially  viscous  type,  with  strain  rates 
written  as  functions  of  stress  only,  and  vice  versa.  Other  parameters  such  as 
plastic  strain,  plastic  work,  temperature  etc.  are  treated  as  parameters.  These 
can  easily  be  included  or  changed  as  may  be  necessary  for  particular  estima- 
tion purposes.  The  use  of  essentially  viscous  constitutive  relations  in  the 
sense  defined,  implies  that  plastic  rate  sensitivity  is  the  primary  material  be- 
havior. These  relations  include  perfectly  plastic  behavior  as  a special  case. 

The  situation  is  exactly  analogous  to  the  use  of  perfectly  plastic  behavior  as 
the  basis  for  plastic  limit  analysis  in  conventional  structural  theory;  without 
this  idealisation  the  limit  theorems  would  not  exist,  and  this  powerful  tool  for 
analysis  and  design  would  not  be  available. 

A further  step  in  the  application  of  the  basic  theorem  is  the  introduction 
of  homogeneous  relations  of  viscous  type.  Strain  rate  test  data  are  conveniently 
plotted  as  curves  of  stress  at  fixed  strain  levels  as  function  of  strain  rate 
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on  a logarithmic  scale.  For  the  three  structural  metals  mentioned,  and  some 
others,  the  curves  are  concave  upward.  The  experimental  curves  then  can  be 
reoresented  verv  closely  bv  equations  of  the  form 

• 1/n  a i 

a = a + be  or  — = 1 + (— ) (la,b) 

where  a > 0 , e > 0 ore  a uniaxial  stress  and  coniu^ate  strain  rate,  respec- 
tively, and  a,  b,  n (or  equivalently  a , £ , n)  are  constants  tor  a c^iven 

o o 

strain  level.  More  convenient  for  manv  purposes  are  honoeeneous  forrs.  (?or 
example,  the  stress  path  for  maximum  complementary  work  takes  a particularly 
simple  form  when  the  stress-strain  rate  relations  are  homogeneous . ) If  the 
constants  in  Eqs.(l)  are  determined  from  experiments,  a homogeneous  relation 
can  be  derived  from  them  in  the  form  [id] 


where  n — „ - 

n - - fll/n 

and  3 = e/eo  » c beinr.  the  strain  rate  at  which  the  expressions  (1)  and  (2) 
are  "matched',  in  the  sense  that  the  two  curves  touch  with  common  tangent  at 
this  point.  The  homogeneous  form  (2)  alwavs  lies  below  the  inhomoqeneous  form 
(1);  hence  the  replacement  of  (1)  bv  (2)  leads  to  a consistent  error  in  deflec- 
tion estimate  of  conservative  type.  Comparison  of  solutions  for  one  and  two 
decree  of  freedom  problems,  and  for  the  few  continuous  structures  for  which 
solutions  are  available,  have  shown  extremely  small  errors  when  the  matching 
strain  rate  I is  representative  of  the  highest  strain  rates  occurring  in  the 
response.  The  determination  of  the  inhomogeneous  form  (1)  from  test  data  is 


1/n' 


(f-  ) 


(2) 


j i 

'o 


1 + 3 

a “ u 1/n' 
o 3 


,1/n 


( 3a ,b) 
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discussed  in  [2],  together  with  comparisons  between  final  deflections  and  re- 
sponse times  from  the  two  types  of  representation.  3oth  are  readily  written  for 
general  combined  stress  states  and  in  terns  of  appropriate  generalized  stress 
and  conjugate  strain  rates.  This  is  done  in  various  forms  in  references  [1-9]. 
Effects  of  strain  rate  history  are,  of  course,  not  considered  in  either  form. 

The  application  of  the  theorem  requires  the  solution  of  a problem  of  static 
loading.  To  obtain  an  upper  bound  on  the  deflection  at  a certain  point  and  in 
a specified  direction,  a force  at  this  point  parallel  to  the  stated  direction  is 
required  to  induce  total  strain  enerspy  at  least  as  larr,e  as  the  initial  kinetic 
enertpy  of  the  dynamic  problem.  This  was  approached  bv  iterative  schemes  in 
[1,  2,6].  In  [6]  the  equations  of  the  bound  calculation  were  written  in  forms 
identical  to  those  of  the  equations  involved  in  the  mode  iteration,  with  dif- 
ferent meanings  of  constant  parameters. 

A difficulty  arises  because  of  the  presence  of  the  response  time  of  the 

dynamic  problem  in  the  basic  inequality  of  the  method.  This  can  be  replaced  by 

an  upper  bound  on  this  time.  Although  both  upper  and  lower  bounds  on  response 

time  can  be  calculated  by  known  methods,  these  are  valid  under  assumptions  of 

"small  deflections"  and  hence  not  strictly  applicable  to  the  present  problems. 

However,  the  deflection  estimate  is  extremely  insensitive  to  the  value  of  re- 

1/n ' 

sponse  time,  which  appears  in  the  form  t^  , where  n'  for  the  structural 
metals  of  interest  here  is  at  least  5 and  may  be  as  larr^e  as  15  or  20. 

Experiments  are  necessary  to  investigate  the  practical  validity  of  com- 
puted deflection  bounds;  the  calculation  is  not  exact,  but  the  sitpis  of  errors 
due  to  the  idealizations  and  simplifications  made  can  be  assessed.  Experiments 
were  therefore  made  on  small  frames  of  mild  steel  and  titanium,  including 


material  behavior  tests  at  moderately  high  rates.  These  experiments  are  de- 
scribed in  detail  in  reference  [71.  The  tests  were  fairlv  comnrehensive . For 
both  materials  two  tvpes  of  frame  were  tested,  corresponding  to  impulsive  load- 
in  p;,  concentrated  or  distributed  over  the  beam  member.  For  each  material  and 
frame  loading  tvpe,  impulses  were  aoolied  of  magnitudes  over  a wide  range,  so 
that  the  maximum  permanent  deflections  ranged  from  elastic  magnitudes  to  about 
a third  of  the  span  ( 40  or  more  thicknesses). 

The  calculated  deflection  bounds  were  found  in  all  cases  to  lie  above  the 
observed  final  deflections  of  the  tests.  The  major  idealizations  and  simplifi- 
cations were  deliberately  of  a consistently  conservative  nature  (i.e.  such  as 
to  increase  the  deflection  bound).  The  bound  magnitudes  were  found  to  exceed 
the  test  values  by  factors  as  large  as  50  to  75  nercent.  In  view  of  the  sim- 
plicity and  cheapness  of  the  calculation  compared  to  a full  numerical  solution, 
over-estimates  of  this  extent  would  be  entirelv'  acceptable  for  most  practical 
purposes . 

Although  the  calculation  required  is  certainly  "relatively  simple",  in 
its  form  as  illustrated  in  references  [l,  2,  6]  the  basic  simplicity  is  to  some 
extent  obscured  bv  manipulations . Work  is  now  proceeding  (under  other  sponsor- 
ship) to  make  the  calculation  more  direct  and  systematic,  in  order  to  enhance 
the  practical  utility  of  the  method. 





Mode  Approximations 
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The  so-called  "mode  approximation  technique  1 was  proposed  by  Martin  and 
Symonds  [11]  for  impulsively  loaded  structures.  The  basic  theory  was  es- 
tablished for  rigid-perfect lv  plastic  behavior  and  for  small  deflections.  It 
was  later  shown  by  Lee  and  Martin  [12,  13]  and  Lee  [14,  15]  that  structures  of 
other  material  behaviors,  in  particular  viscoplastic,  could  be  treated  in  the 
same  scheme.  The  concept  of  "instantaneous  mode  solution"  was  shown  to  be 
available  for  viscoplastic  behavior  which  because  of  the  constitutive  equations 
are  nonhonoqeneous , does  not  permit  mode  form  solutions  to  persist.  The  pro- 
cedure was  limited  to  small  deflections,  and  even  for  these  was  rather  lengthv. 

Initial  contributions  under  the  present  research  program  [2]  showed  a wav 
of  applying  the  basic  concepts  of  the  mode  approximation  method  to  impulsively 
loaded  structures  reaching  large  deflections,  including  the  important  cases 
where  the  response  changes  qualitatively  at  finite  deflections,  and  with  a pos- 
sibility of  extension  to  more  general  pulse  loading  and  material  behavior.  At 
the  same  time  the  handling  of  strongly  rate  sensitive  plastic  behavior  v;as 
shortened  by  the  use  of  the  homogeneous  "matched  viscous"  representation  [10]. 
As  already  pointed  out  above,  this  simplification  is  not  only  realistic  but 
consistently  conservative. 

The  proposed  new  approach  makes  use  of  special  "natural  motions"  in 
separated-variable  (mode)  form.  These  satisfy  the  full  system  of  field  equa- 
tions (dynamics,  compatibility,  constitutive  equations)  but  do  not  in  general 
agree  with  specified  initial  velocity  values  due  to  impulsive  loading.  The 
existence  of  persistent  ("permanent")  solutions  of  this  form  requires  certain 
conditions  to  be  met,  in  particular  linearity  of  the  equations  of  dynamics  and 
compatibility  and  homogeneity  of  the  constitutive  eouations.  When  these  con- 
ditions are  not  satisfied,  as  in  large  deflection  problems,  instantaneous  mode 
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form  solutions  can  still  be  identified  for  a particular  deflection  field  and 
measure  of  velocity  -agnitude.  The  amplitude  of  the  initial  mode  form  solution 
is  chosen,  as  in  the  original  small-defiection  version,  so  as  to  minimize  the 
initial  value  A°  of  the  functional 


A[u^(x,t  ) 


U.(x,t)]  = x-  / o(u. 
1 i v 1 


u.)(u.  - u.)dV 


(4) 


where  u^ , (i  = 1,2,3)  is  the  actual  velocity  field  and  u.  is  that  of  a 
mode  form  solution;  o is  mass  density,  and  the  integral  covers  the  structure. 
The  convergence  property  (dA/dt  $ 0)  can  no  longer  be  proved  for  large  deflec- 
tions, but  it  may  still  be  postulated  that  the  "best"  initial  mode  field  re- 
mains the  one  which  minimizes  the  initial  mean  square  difference  between  the 
specified  and  the  mode  field,  with  mass  densitv  as  a weighting  function.  This 
provides  a simple  and  advantageous  starting  condition,  which  is  generally  better 
than  others,  for  example  matching  initial  kinetic  energy. 

The  new  aporoacn  was  applied  first  to  a simple  discrete  structural  model 
[2]  and  then  to  a circular  plate  [3,  3]  and  rectangular  frame  [6].  The  main 
problem  in  these  applications  is  the  determination  of  the  mode  velocity  field 
from  the  field  equations,  together  with  the  associated  acceleration  magnitude. 
This  nonlinear  eigen  problem  was  solved  by  iterative  schemes.  In  [3]  this  was 
done  making  use  of  a finite  element  format. 

Comparisons  with  experiments  are  essential  to  test  the  validity  of  such 
an  approach,  which  makes  use  of  idealizations  and  simolifi cat ions  and  which  in- 
volves also  a definite  "intrinsic"  error  (in  determining  the  starting  velocity 
amplitude  by  the  "minimum  A°"  technique.)  Hence  tests  were  carried  out  on 
circular  plates  and  frames,  designed  to  investigate  both  tyres  of  error.  The 
tests  on  frames  [7]  have  already  been  mentioned.  Those  on  plates  were  also 
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quite  comprehensive.  Two  rate  sensitive  metals  were  used  (mild  steel  and  ti- 
tanium), and  three  loading  conditions,  differing  in  the  degi>ee  of  concentration 
of  the  impulsive  pressures.  The  specimens  were  "full/  clamped"  circular  plates. 
For  each  of  the  six  cases  a range  of  total  impulse  magnitudes  was  applied  so 
that  final  permanent  displacements  ranged  from  elastic  magnitudes  to  about 
eight  plate  thicknesses.  Full  details  of  the  tests  are  given  in  [9]>  including 
the  determination  of  constants  for  strain  rate  sensitive  behavior  and  tech- 
niques for  measurement  of  impulse  and  time  history  of  displacement.  The  ap- 
proximate determination  of  response  by  a sequence  of  instantaneous  modes  is 
described  in  [9].  Here  the  calculation  is  somewhat  improved  in  efficiency  over 
the  earlier  version  [3],  by  making  use  of  "master  response  curves"  for  struc- 
tures of  given  geometry  and  class  of  material  behavior.  An  iterative  scheme 
to  determine  instantaneous  mode  velocity  fields  was  still  employed,  but  this 
was  cast  in  finite  element  format  for  convenience.  Convergence  difficulties 
were  encountered  at  large  deflections,  but  from  empirical  observations  these 
are  believed  to  be  of  no  importance  as  far  as  final  deflections  are  concerned. 

The  comparison  of  test  results  (final  maximum  deflections,  deflection  pro- 
files, and  response  times)  with  predictions  of  the  mode  approximation  technique 
are  discussed  in  [9].  Here  emphasis  is  put  on  identifying  the  main  sources  of 
error.  The  agreement  between  test  deflection  magnitudes  and  those  predicted 
is  in  general  quite  good,  although  not  so  close  as  found  in  the  frame  tests  (see 
C6,  7]).  This  agreement  is  less  important  for  present  purposes  than  understand- 
ing why  good  agreement  should,  or  should  not,  be  expected.  In  both  the  frame 
and  the  plate  test3,  some  anomalous  results  remain  to  be  explained,  as  discussed 
in  [6]  and  [8].  Although  the  "clamped  plate"  test  is  a prototype  of  the  type 


of  structure  whose  response  changes  qualitatively  as  large  deflections  are 
reached,  the  full  clamping  condition  is  experimentally  difficult  to  achieve. 

It  should  be  emphasized  that  the  "mode  technique'1  is  not  merely  a device 
for  treating  dynamic  plastic  structural  response  by  one  degree  of  freedom  models 
and  ad  hoc  short-cuts.  The  engineering  literature  is  full  of  approaches  of 
this  nature;  by  far  the  most  rational  and  general  such  approach  is  that  of 
Kaliszky  [16,  17].  The  concepts  of  mode  form  response  for  structures  of  certain 
classes  of  inelastic  behavior  are  as  fundamental  in  dynamic  Dlastic  theory  as 
are  the  concepts  of  limit  load  and  plastic  collapse  in  static  structural  analysis. 
Close  links  exist  between  limit  analysis  and  dynamic  plastic  mode  response. 

Some  of  these  were  suggested  by  Lee  [14]  and  Martin  [18].  Our  research  program 
included  contributions  to  understanding  these  links,  both  from  a direct  point 
of  view  illustrated  by  structural  examples  [5],  and  from  that  of  general  varia- 
tional-extremal properties  which  characterize  mode  form  responses  [4].  The 
latter  discussion  was  motivated  by  the  hope  of  improving  practical  computational 
schemes  by  the  availability  of  well  understood  variational-extremal  theorems 
of  both  kinematic  and  dynamic  type. 
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Summary 

The  research  completed  and  described  in  the  nine  reports  or  publications 
of  References  1-9,  has  made  a contribution  to  the  development  of  a branch  of 
basic  theory  underlying  response  of  structures  to  severe  pulse  or  impact  load- 
ing [2,  4,  5].  This  theory  has  direct  application  to  rigorously  based  general 
methods  for  obtaining  estimates  of  the  final  maximum  deflections  and  response 
times,  as  exemplified  by  [1,  2,  3,  6,  8].  Rather  comprehensive  series  of  ex- 
periments were  carried  out  [7,  9]  on  structures  of  tvpes  chosen  to  exhibit  dif- 
ferent large  deflection  behaviors  of  particular  interest.  The  generally  close 
agreement  between  predicted  and  observed  deflections  and  times  was  critically 
discussed  in  [6,  8]  with  emphasis  on  the  sources  of  error. 

A considerable  increase  in  understanding  of  some  fundamental  features  of 
dynamic  plastic  response,  as  contrasted  with  facility  in  producing  sheets  of 
computational  print-outs,  may  be  said  to  have  been  achieved. 

A number  of  important  questions  remain  unsolved,  concerned  particularly 
with  implementing  the  approximation  technioues  in  a more  standard  and  systematic 
fashion.  These  are  now  being  pursued  under  new  auspices  with  more  adeouate 
funding. 
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